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Synopsis 


Laminated composite plates are fast replacing metal alloys in most light 
transport vehicles and various other industries. Increasingly, many aerospace and high- 
speed rail components are being fabricated with composites. 

The material properties of such composite display considerable scatter because of 
the uncertainty involved at many levels-properties of its constituents, fabrication and 
manufacturing processes, geometrical parameters of laminates, fiber orientations, volume 
fraction, curing temperature and pressure, voids, impurities, curing time etc. This leads to 
variation in stiffness coefficients of the laminates and thus results in uncertainty in the 
response behavior. Besides this the structure may be subjected to variations in input 
loading. The system can best be modelled as having random material properties with 
random excitation. In some applications the structure undergoes large deformations, 
which demands the investigation of their nonlinear behavior. 

The present work is an attempt * to incorporate nonlinear effects in the random 
environment with the material properties and external loading as random. The second 
order statistics of static deflection, natural frequencies and forced vibration response of 
composite plate with random material properties and random excitation is evaluated. The 
basic formulation uses classical laminate theory and Von-Karman non-linear strain - 
displacement relationship. The system equation for the static case is obtained by using 
Rayleigh-Ritz method whereas the dynamic formulation uses Hamilton’s principle. The 
numerical results for mean and SD for the responses have been presented. 
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CHAPTER 1 
INTRODUCTION 

1.1 GENERAL INTRODUCTION 

Composites are known to have specific advantages over metals as these can be 
tailored to have directional properties required for the loading. This has lead to their 
extensive use in diverse engineering applications. Composite has a large number of 
parameters involved with its production and fabrication process. Change in factors such 
as fiber orientation, laminae thickness, volume fraction, curing temperature and pressure, 
voids, impurities, curing time etc. induce variations in the lamina properties. This leads to 
variation in stiffness coefficients of the laminates. It is practically not possible to 
completely control variation in all these parameters and thus dispersions in the material 
and geometric properties are inherent. In use, the uncertainty in the system parameters 
result in uncertainty in the response behavior. Further, the external loading (both static 
and dynamic) may also be random in nature. In many applications, like aerospace 
engineering, it is important for the designer to have an accurate evaluation of the 
structural response. For such sensitive applications, enhanced accuracy in response 
evaluation is possible by modeling the problem with random material properties and 
random external loading. Some of these structures are very often subjected to severe 
loadings that result in large response and consequently demand the investigation of their 
nonlinear behaviour. 

1.2 LITERATURE SURVEY 

Any structural analysis problem is characterised by the following three basic 
aspects: material of the structure, geometry of the structure and type of loading. In real 
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life problems, all the three aspects are random in nature. Significant volume of literature 
is available for external loading as random with material properties and geometry as 
deterministic. Nigam and Narayanan [1] present various classes of problems in this area. 

Some published literature is available for analysis of conventional structure with 
random material properties. Ibrahim [2] has reviewed topics pertaining to structural 
dynamics with parameter uncertainties. Nakagiri et al. [3] have studied simply supported 
(SS) Graphite/Epoxy plates with stochastic finite element method (SFEM) taking fiber 
orientation, layer thickness and layer numbers as random variables and found that the 
overall stiffness of fiber reinforced composite (FRP) laminated plates is largely 
dependent on the fiber orientation. Fukunaga et al. [4] have investigated the effects of 
scatter in the lamina strengths, relative fiber volume fraction, and laminate stacking 
sequence on the ultimate strength of the hybrid laminates. Raj et al. [5] have analysed 
rectangular plates with and without cut-outs using higher order shear deformation theory 
by using combination of FEM and Monte Carlo Simulation (MCS). 

Leissa and Martin [6] have analysed the vibration and buckling of rectangular 
composite plates, and have established that variation in fiber spacing tends to increase the 
buckling load by 38% and the fundamental frequency by 21%. Shinizuka and Astill [7] 
have employed a numerical technique to obtain statistical properties of eigen values of 
spring supported columns with the spring support and axial loading along with material 
and geometric properties as random. The method has been used to investigate the 
accuracy of perturbation approach for calculation of vibration and buckling modes. Salim 
et al. [9-12] have studied the statistical response of plates considering material properties 
as independent random variables (RV). The second order statistics for static deflection. 
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natural frequency and buckling load of rectangular plates have been studied using a first 
order perturbation technique (FOPT). Free vibration response has been obtained by 
Vaicatis [13] for beams with mass and flexural rigidity as random variables. The 
nonlinear forced vibration of circular and rectangular plate with various boundary 
conditions has been studied [15-18] by applying the Galerkin or Rayleigh Ritz methods. 
Singh et al. [19] has studied the nonlinear forced vibration of antisymmetric rectangular 
cross-ply plates by using Hamilton principle. 

Chen and Soroka [20] have studied the response of a multi degree of freedom 
system with random properties to deterministic excitations. The system equations have 
been solved by perturbation technique. The second order statistics of the system response 
have been investigated with variation in the system property statistics. Yadav and Verma 
[21-22] have studied the buckling and free vibration response of thin cylindrical shell 
using classical laminate theory (CLT) and have employed the FOPT for obtaining the 
second order statistics of buckling loads and natural frequencies. Gorman [23] has 
presented free vibration analysis of thin rectangular plates with variable edge supports 
using the method of superimposition. Singh et al. [24-25] have studied the initial 
buckling and natural frequency of cylindrical panel and composite plate with random 
material properties and have obtained the second order statistics of response. 

PRESENT WORK 

Many researchers have studied the uncertainties in material strength and stiffness of 
composites. A comprehensive survey of the literature indicates that a significant volume 
of publication exists for response of composite beam and plates with linear strain- 
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displacement relation. However, the response of composite plate in random environment 
with nonlinear strain-displacement has not received attention to the researchers. In many 
applications, the deflections are large and may lead to nonlinear behaviour. The present 
work is an attempt to incorporate the nonlinear effects in the random environment. The 
second order statistics of static deflection, natural frequencies of free vibration and forced 
vibration response of composite plate with random material properties and random 
excitation is evaluated. The basic formulation uses CLT and Von-Karman non-linear 
strain displacement relationship. The system equation for the static case is obtained by 
using the Rayleigh-Ritz method whereas Hamilton’s principle has been employed for the 
dynamic problems. Lamina material properties are modeled as basic random variables 
(RVs) and are assumed to be uncorrelated for the generation of the results. External 
excitation is assumed to be random in time. The numerical results for mean and SD for 
the responses for known statistics of material properties and loading are obtained. Effects 
of parameters like side-to-thickness ratio and aspect ratio for symmetric and anti- 
symmetric cross-ply laminates has been examined. 
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CHAPTER 2 


FORMULATION 

2.1 INTRODUCTION 

In this chapter expressions are developed for the response statistics of composite 
laminated plate with random material properties having simply supported edges. Three 
types of responses are discussed in this study- (i) static deflection, (ii) free vibration and 
(iii) force vibration. As mentioned before, the classical laminate theory and Von Karman 
non-linear strain displacement relation have been used for the basic formulation. The 
system equation for the static case is developed by using Rayleigh Ritz method whereas 
for the dynamic cases Hamilton principle has been employed. First order Perturbation 
technique is used to analyse the problems and the second order response statistics has 
been obtain. 

2.2 GENERAL FORMULATION 

Consider a rectangular composite plate of constant total thickness l h\ composed 
of thin orthotropic layers bonded together. The origin of a Cartesian coordinate system, 
as shown in Figure 2.1, is located in the central plane at the left comer with x and y axes 
along the middle plane and the z-axis normal to this plane. 

If Qij are the reduced stiffness matrix terms, the simplified relations between basic 
lamina properties and reduced stiffness matrix terms are given as [26]: 

Qu—E\\K\-Vi\ H 2 ) 

Qll-Elllil-VlX V\j ) 
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Q\2-V\tEi«iI(\-Vi\ V12) 

Qi\-Q\i and Q66~Gn 


( 2 . 1 ) 


V 21 can be expressed in terms of the other material properties as: 
v lx =v n E 12 IE n 

The transformed reduced stiffness matrix terms Qy are related to the reduced stiffness 
matrix terms Qy through the following equations [26]: 

Q 1 \=Q\\cos a 9 +Qi 2 sin 9 +2(Qn + 2Q66)sin 2 9 cos 2 6 

Q 2 i = Q\ 1 sin 4 6 +Q 22 COS 4 9 ^' 2 .{Q\ 2 Ar ' 2 Q(,^)siti' 9 cos 2 9 

Q \i=Qn{cos 4 9 +sin 9 )+(Q\ \+Q22-4Q 6 6)sin 2 9 cos 2 9 (2.2) 

Q 66 = Q66(cos 4 9 +sin 4 9 )+(<2i x+Qu-'ZQn-'ZQ^sin 2 9 cos 2 9 

Q \6 = (Qn-Qn-2Q66)cos 3 Osin 9 -(Qi2-Q\2-2Q( > ( 1 )sin 3 9 cos 9 

Q 26 = (Q\\-Q\ 2 - 2 Q 66 }cos 9 sin 3 9 -{Qr2-Qn-2Q(^)sin 9 cos 3 9 
The elements of the extensional stiffness matrix Ay, coupling matrix By and bending 
stiffness matrix Dy are given by [26]: 

A ij= Z (Q ij)k{hk -h-\) 

k = 1 

5ij=(l/2) X (Qij)k( h k (2-3) 

k = 1 

Aj=(l/3) Z (QiHK -hL) ij = 1,2 and 6 

*=1 

where ‘A*’ is the thickness of the tf h lamina in the given laminate. 
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The stress and moment resultants per unit length of the plate in terms of stiffness matrices 
are defined as [26]: 




[ A S 

V 

V 

Mi. 


I 


_ K J_ 


(ij=l,2,6) 

where e . and Kj are the mid plane strains and curvature strains respectively. 
The strain energy for the plate can be written as [3 1]: 

l hi 2 

u =~ J K + + & 

2- -hi 2 A 

The stress and moment resultants can also be defined as [31]: 

(i=x.y) 


(2.4) 


(2.5) 


A/2 


(N it Mi)= | a i (1, z)dz 

-A/2 


( 2 . 6 ) 


Neglecting transverse shear effect, i.e. (y xz ,y yx = 0) and by substituting equations (2.4) 
and (2.6) in (2.5) the expression for the strain energy of the plate becomes: 

\ a b 

U = —\\{A\\s x +2A\2 S x s +A 22 s +2A\ 6 s x / +A 26 S/ +A66/„ + 


2 00 


2B\ 1 k x s x +B n (K y e x + 2K x £ y )+2B 22 K y £ y +2B l6 (K xy s x + Kj xy ) + 


2B 16 {K xy s y + >c y Y xy ) + 2B 66 y xy K xy +D n e x 2 +2D n K x K y +D 12 £ y 2 + 


2D 2 (,rc xy K y +2Di 6 K xy K x +D 66 K xy 2 }dxdy 
Assuming that the deflection may be large •* 
strain-displacements are used [31]: 

Sv ! ,dw. 2 


(2.7) 

, the following Von Karman 


du | 1 
dx 2 dx 


»,-* + *<** %-■=:+-+ 


du + dv + dw dw 
dy dx dx dy 


( 2 . 8 ) 


where e x , s and y v are mid plane strains and u, v and w are mid plane displacement in 
x, y and z directions. 
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The mid plane curvature-displacement relationship can be put as: 


K 


X 


d 2 w _ d 2 w _ 0 d 2 w 
US’ Ky ~~~dy r ' Kxy ~ tedy 


(2-9) 


2.3 BOUNDARY CONDITION 

The boundary condition for a^plate with all sides simply supported with edges free 
to move in their respective normal directioitfare given by the following set of equations: 
Along x=0 and x-a for all y, 

d 2 w d 2 w 

w=0: 

ox oy 

v=0: N x -A\\ — +A\2~=0 
dx dy 

Along y =0 and y-b for all x, (2.1 0) 

d 2 w d 2 w 

w=0: M y = -D n { —y )-D 22 ( ~y)=0, 
ox oy ' 

u= 0: N y =A n ^ + A 22 ^=0 
ox oy 

2.4 STATIC DEFLECTION ANALYSIS 
2.4.1 System Equations 

Consider a specially orthotropic rectangular laminated composite plate subjected 
to a distributed sinusoidal transverse loading given by: 

Q(x,y)= Q 0 Sin — Sin ^ (2.11) 

a b 

where a and b are the plate dimensions in x and y directions respectively, Q 0 is a constant 
giving the loading amplitude and m and n are integers. This loading form is generic and 
can be used to simulate different types of transverse loadings on the plate. 
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In case of rectangular antisymmetric crossply all the coupling elements of [A], 

[B] and [D] matrix identically go to zero i.e. A i 6 =A 26 =B i (,=B 2 6=D 1 6 =JD 2 6=0 . From 

equation (2.7) the strain energy for the plate can be written as: 

1 a b 2 2 2 

U=- )\{A u e x +2Ai2£ x s y +AnS y +A 66 y xy +2B n k x s x +B n (K y £ x +2/c x s y ) 


+ 2B 2 2K y £ y +2B 66 y xy K xy +D u s x 2 +2Dx 2 K x K y + D 22 £ y 2 + D 66 K xy 2 }dxdy (2.12) 
For a sinusoidal static load defined in equation (2.11), the external work done can be 
represented by: 

W= J | Q(x,y) w(x,y) dxdy (2.13) 

oo 

The following set of admissible functions can be used to satisfy the defined boundary 
conditions: 

^ mi vc njty mux _ nny 

u=uqCos Sin , v=vo Sin Cos — — , 

a b a b 

mux nity 

w=woSin Sin — — (2.14) 

a b 

where u 0 , v 0 , and wo are constants indicating the displacement amplitudes. 

Substituting equations (2.8), (2.9) and (2.14) in equations (2.12) and (2.13) and 

integrating gives the total potential energy of the system as: 


K=U-W (2.15) 

Minimizing the total potential with respect to mo, vo and vv 0 results in the following three 
simultaneous nonlinear algebraic equations: 

Txuo+Tyvo+Tiwo+Ttwo 2 ^ (2.16) 

T 2 uo+T5Vo+TeWo+TjWo 2 =0 (2.17) 

r 3 uo+T 6 vo+T & wo+2T 4 uoWo+2T 1 voWo+T 9 wo 2 +TxoW ~Qo (2. 1 B) 
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where the coefficients T\, T 2 ,—,T\o depend on the plate geometry, material properties and 
the mode shapes. Their expressions are: 

T\={mn /a) 2 A \ i +(« n / Vf A T 2 =(m n /a){n n fb){A 1 2 +A 6 q) 

Ti= -{mn/dfB\\, 

T A =(-4/9mn n 2 )S mn [(m n /a) 3 A \ \+(m n /a)(n n /b) 2 {A\ 2 -A^)} 

T s ={m tc /a) 2 A 66 +(n jt /b) 2 A 22 , T 6 = -{n n /b) 3 B 2 2 (2. 1 9) 

Ti =(-4/9mn n 2 )S mn [(n n /b) 3 A 2 2 +(n n /b){m n /a) 2 (A 1 2 -A 6 a)] 

Tz=(m tt /a) 4 D\ 1 +2(m n /a)\n n /bf{D\ 2+2D 66 )+(« n /b) 4 D 22 
T 9 =(4/3mn n 2 )S m „[(m n /a) A B\ 1 +(n n /b) 4 B 22 

T\o=(9/32)[(m n /a) 4 A 1 1 +(n n /b) 4 A 22 ] +(1/1 6 ){m n /a) 2 {n n /b) 2 (A 1 2 +2A 66 ) 
where iS m „=(l-(-l) m ) (l-(-l)") 

It is possible to obtain u 0 and v 0 in terms of w 0 from equations (2.16) and (2.17). 
Thereafter, substituting in equation (2.18), the resulting equation in w 0 becomes: 

«o=[ {(T 2 T 6 -T 3 T 5 ) w 0 +(T 2 T 7 -T 4 T 5 ) w] }/(T,T 5 -T 2 2 )] (2.20) 

vo=[{(T 2 T 3 -TtT 6 ) w 0 +(T 2 T 4 -T,T 7 ) w 2 }/(T,T 5-T 2 2 )] 

and L x w 0 +L 2 w 2 +Z 3 w 0 3 =g 0 (2.2 1 ) 

Where L\,L 2 and L 3 are as follows: 

Li=[t s + {(2r 2 r 3 r 6 -r 3 2 r 5 -r 6 2 r,y(r 1 t s -t 2 )}} 

L 2 =[T 9 + {(3 T 2 T 2 T 2 +3 T 2 T a T 6 -3 T 3 T a T 5 -3 T, T 6 T 2 y(T x T 5 -T 2 2 )}] (2.22) 

Li=[T l0 + {(4T 2 T a T 7 -2T 5 T a 2 -T x tMTx T 5 -T 2 2 )}] 

2.4.2 Mean and variance of static response - perturbation technique 

First order perturbation approach has been adopted for obtaining the deflection 
statistics with randomness in material properties and loading on a flat plate. Here it is 
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assumed that all the material properties and loading components are independent of each 
other. Also it is assumed that the dispersion of each random variable about its mean value 
is small, which is true in most engineering applications. 

Any random variable may be split up as the sum of the mean variable and zero 
mean random part with generality as: 

R=R d +R\ (2.23) 

where superscript ‘d' denotes the mean value, which is deterministic, and V denotes the 
superimposed zero mean random component. 

Using Taylor series expansion the random part of the dependent variables can be 
expressed in terms of the independent variables. Consider L\ as the system operator, w as 
the displacement response of the system, Q to denote the force applied to the system and 
bi as the basic material properties. The primary variables are assumed to be the basic 
material properties b[ and the elements of the applied load < 7 *. As assumed, the random 
part in the primary variables is small in magnitude compared to their mean values, the 
second and higher order terms are neglected and the expression may be put as: 

L r i ='Z(dL d /dbf)b{ 

1 

<=Udw d 0 /dbf)b;+ Udw d 0 / dqf )q[ (2.24) 

/ k 

Since dispersion about mean is small and only first order terms are retained, derivatives 
of the random process have been approximated as the derivatives of their deterministic 
values. 

The basic random variables 6 / have been sequenced as follows: 

b\=E\\, b 2 =E 22 , by=v n and bf=Gn (2.25) 
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Using Taylor series expansion, keeping only one term in series and neglecting small 
quantities one can write the random part of reduced stiffness matrix and transformed 
reduced stiffness matrix as: 

Qs=i(s(Q^yabf)b; and Q;-tmQiyaf)b; (2.26) 

7= 1 7=1 

The random part of extensional, bending-extension coupling and bending stiffness 
matrices can also be expressed similarly by using Taylor series expansion keeping only 
one term in the series and neglecting all higher order terms as: 

4 =t(d( Afj ydbf ) b \ , b: = ± (3 ( By )/dbf ) b \ , 

/= 1 1 

and Dy = t(d(Dfj)/dbf)b; (2.27) 

7=1 

where l p' denotes the number of basic RV’s. 

The partial derivative of the reduced stiffness matrix with respect to each material 
property can be written as: 

(i) Taking E\ 1 as a random parameter 

dQn / dE^ =E u (Eu- 2E 22 vf 2 )/(E u -2E 22 v 2 2 ) 2 
8Q d 22 /dE^-E 2 2 v\ 2 /(E\ 1 - 2 E 22 vf 2 ) 2 
dQn / dE^ = - E\ 2 v\ 2 l{En-2Env\ 2 f 
BQ d JdE d n = 0 

(ii) Taking E 22 as random parameter 

QQtJdE d 22 =E 2 n v 2 2 /(E n -2E22V 2 2 f 
dQ d 2 / dE d n = E 2 n l(E \ , -2E 2 2 v 2 n f 
dQ d 2 / dE d 2 =E 2 , v n /(Ei ,-2£ 22 v 2 2 f 
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dQ d JdE d 22 =0 

(iii) Taking v i2 as random parameter (2.28) 

dQf t /dvf 2 =2 £ 2 E 22 v n !{E x l -2E 22 vf 1 f 
dQ d n 1 d 4 =E\ 1 Eh 4 /(£, , -2£ 22 v* ) 2 

30,' / 5 4 = En Ei 2 +E\ ,£ 22 4 /(£, , -2£ 22 v* ) 2 

dQ d Jdvh=Q 

(iv) Taking Gi 2 as random parameter 

8Q d n ldG d n = 0 
dQ d JdG d n =Q 
dQ d 2 /8Gf 2 =0 
8Q d JdG d 2 = 1 

The partial derivatives of the transformed reduced stiffness matrix a ( Q d )/dbf can be 
found from equation (2.28) by using equation (2.2). Also the partial derivatives of the 
stiffness matrix elements A d - , B?j , D d with respect to ‘ bf ’ are required for solving 
equation (2.27). These can be expressed as follows: 

a ( A d ydbf = i(d(Q d ) k /dbf )ih k - h kA ) 

J k =\ j 

a ( B d )/dbf =1 /2 i (a ( Q/ ydbf )( hi - hU ) (2.29) 

d(D d ydbf =1/3 i{d(Q d ydbt){hl -hi,) 

k = 1 

Substituting equations (2.23) and (2.24) in equation (2.21) and collecting terms of 
different orders, we get zero order and first order equation as follows: 
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Zero order equation: 

(If w d 0 +L d 2 (w d 0 ) 2 +Li(w d 0 ) 3 )=Qt 
First order equation: 


(2.30) 


( h\ +2 L{ w d +3 L\ ( w d 0 ) 2 )[ E (dwf / dbf )b;+U dwf / dq d k )q[]+ w d 0 E( dL d ! dbf ) bf + 

/ it / 

{<fY{dL d 2 l dbf ) b\ +( < ) 3 1( 5Zj/ 56/ ) Z>; = 0 O ' (2.31) 

/ ; 

Equation (2.30) gives the mean response of the system whereas equation (2.31) 
can be used to find the response variance of the system. In equation (2.31) dwf / dbf and 
dwf / dqf are the unknowns. 

To find dwf /dbf , multiplying equation (2.31) by b u r and taking expectation we 
get the following under the assumption that the primary material variables and the 
loading components are independent random processes. 

{Lf+2L d 2 wf+3Lf(wf) 2 }(dwf/dbf)+(dL d i /dbf)wf+(dL d 2 /dbf)(wf) 2 + 


(dL d 3 /dbfXw d 0 f= 0 (u=l, 2,..., 4) (2.32) 

To find dwf / dqf , multiplying equation (2.31) by qj and taking expectation we 
get the following with the assumption that b? s and qic’s are independent of each other. 

(Lf+2L d 2 wf +3 Lf ( wf ) 2 )( dwf / dqf )=1 («= 1) (2.33) 

Thereafter, the total response can be evaluated as follows: 

™o=wf+[Z(.dwf IdbfWt+YiHldqf)^} (2.34) 

7 k 

The variance of the static response takes the form: 

vaiw-EUXfav; iaq)b;+j.(pw d jaq d t )q[} 2 } (2.35) 

7 k 


15 



The SD is obtained as the square root of the variance. The above approach is 
general in nature, and can be applicable to a variety of problems in analysis of structures 
involving random system parameters and loading. 

2.5 FREE VIBRATION ANALYSIS 
2.5.1 System Equations 

Consider a rectangular composite plate of constant total thickness l h\ composed 
of thin orthotropic layers bonded together and is executing free vibration. 

In case of rectangular antisymmetric crossply all the coupling elements of [A], 
[B] and [D] matrices identically go to zero i.e. A X6 =A 1( =B , 6 =B 2 6=D \ 6 =-D 2 6=0 . From 
equation (2.7) the strain energy for the plate can be written as: 

1 a b 

U=— f\{A \\S X +2An£ x s v +A 2 iS y 2 +A 66/„. 2 + 2B]\rc x e x +B\ 2 (k s x +2k x s ) 

2 oo 

+ 2B 2 2fc y £ y +2B 66 r xy K: xy +D U £ X +2D l2 K x rc y + D 2 i£ 2 +D 66 K xy 2 }dxdy (2.36) 
The kinetic energy of the plate, neglecting in-plane inertia is: 

r=(l/2)JJ( EpA) ^ W) dx dy (2.37) 

00 

Assuming that the amplitude may be large - , Von Karman non-linear strain- 

displacement is used to evaluate the frequency response of the system. 

The following set of admissible functions satisfies the boundary conditions 

defined at equations (2.10): 

, mnx mnx nny 

u=uo(t) Cos Sin — — , v=v 0 ( t)Sin Cos , 

a b a b 

w=w 0 (t)Sin — Sin^L (2.38) 

a b 
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where u, v and w are mid plane displacement in x, y and z directions and «o vo and w 0 are 
maximum displacements at any instant of time. 

Substituting equations (2.8), (2.9) and (2.38) in equations (2.36) and (2.37) and 
integrating gives the total potential energy of the system: 

K=U-T (2.39) 

Use of Hamilton’s principle leads to two algebraic equations and one second-order 
ordinary differential equation in terms of uo, vo and w 0 . 

T\ u 0 +T 2 vo+T 3 wo+T 4 wo 2 =0 (2.40) 

T 2 uo+T 5 vo+T 6 w 0 +T 7 w 0 2 =0 (2.41) 

(IM) % + T 3 uq+ T 6 v 0 + r g w 0 +2 TaUqWq+IT-jVqWq + TgWo 2 + TxqWq =0 (2.42) 

where the coefficients T\, T 2 ,...,T\o depend on the plate geometry, material properties and 
the mode shapes. Their expressions are defined in equation (2.19). 

Substituting uq and v 0 in terms of w 0 , obtained from the first two equations into 
the third equation results in a single second order differential equation containing 
quadratic and cubic nonlinear terms: 

( Z p,h,) W 0 +Li w 0 +L 2 W 2 0 +Li Wq =0 (2.43) 

where L\,L 2 and L 3 are as follows: 

J L I =[r 8 +{(2r 2 r 3 r 6 -r3 2 r 5 -r 6 2 r 1 )/(r 1 7’5-r 2 2 )}] 

i2=[r9+{(3r 2 r377+3r 2 r4r 6 -3r3r47’5-37’ 1 r 6 r 7 y(7'ir 5 -r 2 2 )}] (2.44) 

i3=[rio+{(4r 2 r4ZV-2r 5 r4 2 -7’ir 7 2 )/(r I r5-r 2 2 )}] 

The energy balance equation can be obtained by multiplying equation (2.43) by w and 
integrating with respect to time: 

(Ip, h,) w 0 2 +L X wl +(2/3 )L 2 w l +(1/2 )I 3 w 0 4 =H=const (2.45) 
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When the plate is subjected to free vibration and undergo large amplitude oscillations, the 
constant H in equation (2.45) at wo=w max i.e. vv 0 =0, can be obtained as: 

ff=i, lv ^+(2/3)L 2w ^+ (1/2)1, wL (2.46) 


Substituting this constant H in equation (2.45) yields: 

(ZpA)"o=£i(vvL- w o>(2^^ (2.47) 

In the absence of coefficient Li equation (2.47) at vv 0 =0 has two equal and opposite roots 
±Wmax- Substituting for displacement response w 0 = w max Sin 9 the non-linear time period 
for such a plate can be obtained as: 


2 k 7C ' 2 

T n ,= — =A J 

CO 0 


deJZXpM 


[i,+j4(l+sin ! «)wL] 


(2.48) 


On simplifying, one can rewrite the above equation as: 


„ 4 V2XM) *' r 2 

Inl ?= = J 


dd 


•y/A (1 + c ) o Vl + asin 2 6 


(2.49) 


where c=0.5Li w^/L] and a=c/(l+c) 

Equation (2.49) is in a form of the elliptic integral, which cannot be evaluated in terms of 
elementary functions. An infinite series representation can be generated for the above 
integral by first expanding the integrands in binomial series and then using termwise 
integration: 


t n a 3 , g 4 

nl 22 42 " 82 162 

The non-linear frequency can be expressed as: 

co =2 K /T n l 


•] 


(2.50) 


(2.51) 
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i.e. co =fun(£ u , E 22 , G n , v n and ww) 


(2.52) 


2.5.2 Mean and Variance of Dynamic Response - Perturbation Technique 

Perturbation approach has been adopted for obtaining the second order frequency 
statistics with randomness in material properties of the flat plate. As stated earlier, it is 
assumed that the material properties are independent of each other and the dispersion of 
each random variable about its mean value is small. 

Using Taylor series expansion the random part of the dependent variables is 

expressed in terms of the independent variables as in section 2.4.2. As b r t are assumed to 
be small in magnitude, the second and higher order terms are neglected and the 
expression may be put as: 

z;=Z(5if/5 bf)b; 

i 

a r =Y(do) d /d b^)b; (2.53) 

i 

Since dispersion about mean is small, total random process involved in the 
derivative have been approximated by its deterministic value. Hence from equation (2.52) 
we have 


dco d /dbf = 
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where partial derivative of if , if and a with respect to bf can be expressed in terms of 
E\\, E 22 , V] 2 and Gn as presented for the static analysis. 

The natural frequency can now be expressed as follows: 

co = d +( Z (3 <y d / d bf ) bf (2.55) 

1 

The variance of the natural frequency can be put as: 

Var( to )=E[ { E (3 d / 5 bf ) b\ } 2 ] (2.56) 

1 

The SD is obtained as the square root of the variance. 

2.6 FORCED VIBRATION ANALYSIS 
2.6.1 System Equations 

Consider a rectangular composite plate of constant total thickness ‘ h\ composed 
of thin orthotropic layers bonded together. The plate is subjected to sinusoidal transverse 
dynamic loading given by 

0o(x,y,t)=0o(t) Sin{ — ) Sin{ ^ ) (2.57) 

a b 

where m and n are integers, a and b are the plate dimensions in x and y directions 
respectively and Qo(t) is the time dependent forcing function distribution amplitude. This 
loading form is generic and can be used to simulate different typesof loading on the plate. 

In case of rectangular antisymmetric cross-ply all the coupling elements of [A], 
[B] and [D] matrices identically go to zero i.e. Ai6=yf26 = #i6 = i?26 = A6 = £ ) 26 = 0- From 
equation (2.7) the strain energy for the plate can be written as: 
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I ab 

u= ~ !H a u£ x + 2A n £ x £ y +A 22 £ y +A 66 y xy +2B u k x £ x +B n (>c y £ x + 2tc x £ y ) 

L o o 

+ 2B 22 K y £ y +2B 6 (,y xy K xy +D\\£ x +2D\ 2 K x K y +D 22 £ y 2 +Z>66 k^ 1 }dxdy (2.58) 

For a sinusoidal dynamic load defined in equation (2.57), the external work done can be 
represented by: 

W=JJ Q(x,y,t) w(x,y,t) dx dy (2.59) 

oo 

The kinetic energy of the plate, neglecting in-plane inertia is: 

T=(l/2) J J Cl P ; hi) w 2 (x,y,t)dx dy (2.60) 

oo 

Assuming that the displacements may be large, Von Karman non-linear strain- 
displacement are used to evaluate the system response. 

The following set of admissible functions satisfies the earlier defined boundary 
conditions: 

mm nny mm _ nny 

M=«o(t)Cos Sm — — , v=vo(t)Sm Cos — — , 

a b a b 

w=nv 0 (t)Sin Sin (2.61) 

a b 

where u, v and w are mid plane displacement in x, y and z directions and uq vo and wo are 
maximum displacements at any instant of time. 

The strain energy, work done and kinetic energy can be obtained by substituting 
equations (2.8), (2.9) and (2.61) in equations (2.58)-(2.60) and performing the required 
integration. 

Based on Hamilton’s principle the variation of the Langrangian (U-W -T) provides the 
governing modal equation of motion. 
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( Z Pi h d w 0 +L\ w 0 +L 2 Wq +L 2 m> 1 =Qo(t) 


( 2 . 62 ) 


where L\,Li and L 3 are as follows: 

LHTsH(2T 2 T 3 T6-T 3 2 T 5 -T 6 2 T0/(T^T 5 -T 2 2 )}] 

L 2 =[T 9 +{(3 T 2 T 2 T 7 +3T 2 T 4 T 6 -3T 2 T a T 5 -3Ti r 6 r 7 )/(T, T 5 -T 2 2 ))} (2.63) 

U=[T l o+{(4T 2 T 4 T 1 -2T 5 T 4 2 -T l T 7 2 )/(T ] T 5 -T2 2 )}] 
where T\,T 2 ,...,T\o are coefficients depend on the plate geometry, material properties and 
the mode shape. These expressions are defined in equation (2.19). 

The energy balance equation is obtained by multiplying equation (2.62) by w and 
integrating with respect to time: 

( Z Pi hi) w 2 0 +Li m>1 +(2/3 )L 2 w 3 0 +(1/2)L 3 w 4 0 -Q 0 ( t) w 0 =H=const (2.64) 

When the plate is subjected to a transverse loading Q(x,y,t) and undergo large 
unplitudes, the constant in equation (2.64) at w 0 = , that is, vv 0 =0, can be obtained as 

H=Li wL+ (2/3 )T 2 w^+ (1/2)L 3 w 1 x - Q 0 w mx (2.65) 

Substituting this constant H in equation (2.64), yields: 

( Z P, hd w 0 2 +Li ( HW - w 0 2 )+(2/3)L 2 ( - w 2 )+( 1 /2)L 3 ( - w 0 4 ) 

-(Go w max -go(t) W 0 )=0 (2.66) 

a the absence of coefficient L 2 equation (2.66) at w 0 =0 has two equal and opposite roots 
Wmax- Substituting w 0 = w max Sin 9 and Qo(t) =Qo Sin 6 the nonlinear time period for 
ich a plate can be obtained as: 

__ (2.67) 

0 M,+!t£,(<+™ 2 0) w L- — - 1 

\ 4 W max 

he above equation can be written as: 
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( 2 . 68 ) 


*72 

T n l=4 | 


ddJZiPA) 


o ^/ZJ”-y/[l + c-<f + 6sin 2 0] 

On further simplification, it can be put as: 

Tn| _ 4 VS(M) *j 2 d9_ 

^Ij(l + c-<i) o Vl + asin 2 0 


(2.69) 


where c^OSLiw^/Li, d=Qo/(L\ w max ) and a=c/(l+c-d) 

The above integral is in the term of an elliptic integral, which can not be evaluated in 
.errns of elementary functions. An infinite series representation for above integral is 
generated, as in the previous section, by first expanding the integrands in binomial series 
rnd then using termwise integration. The final expression can be written as 


_ 2^jUp^j g a 2 or 3 a 4 
"'"VAa + c-d) 1 _ 2 2 4 2 8 2 16 2 


(2.70) 


The above equation can also be rearranged to give the unknown w max in terms of the 
mown forcing frequency co . This is a fifth order nonlinear algebraic equation given by: 


A W L +B W Lx +C W Lx +D Wmax+E-0 


(2.71) 


vhere A = 50 L] /64, B-L\L^-Li co 2 Y^ P ; hi, C—(5 0/3 2) (L i L 3 -L 3 P 0 ) , 

D=2 L] -2L\ co 2 YjPi hi and E=2PqLi 
fence the non-linear amplitude can be written as 

i.e. w max = fhn(2?n, E 22 , Gn, V 12 , <0 and Qo) (2.72) 

Equation (2.72) gives the required non-linear amplitude for the specified problem. 
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2.6.2 Mean and variance of forced vibration response - perturbation technique 

Perturbation approach has been adopted for obtaining the deflection or amplitude 
statistics with randomness in material properties and excitation for a flat plat. Here also, it 
is assumed that all the material properties and loading components are independent of 
each other. It is further assumed that the dispersion of each random variable about its 
mean value is small. The random variable may be split up as the sum of mean variable 
and zero mean random part with generality as in equation (2.23) 

Using Taylor series expansion some of the terms of the equation (2.71) can be 
;xpanded as: 

Li-uaLl/dbfw 

i 

<**=£( dw^/dbfW+Zi d»L'dq d k )q r k (2.73) 

l k 


Assuming b\ to be small in magnitude, the second higher order terms are neglected. Here 
• k are RV components of applied loading. Since dispersion about mean is small, total 
andom process involved in the derivative process can be taken as the derivative of the 
eterministic values. 

Using equation (2.72) we get the derivative of w^ x with respect of bf and q k as 
)und in static case. We have: 


dE dA 


dB 


dC 




d _dbf 


db; 


w 

d max 


db, 


w — 

d max 


db\ 


w — 

d max 


dD 


db; 


w 

d max 


5 ^ W max + 3 ' 8W max + 2Cw max + D 


imilarly dw* / dqf can be found as: 


(2.74) 
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(2.75) 




dE__dC_ 2 

ps„d £.d ^nrn 

QQk dq k 


5Aw max + 3^L +2Cw |nJtt +Z) 
where partial derivative of A, B, C, D and E with respect to bf and q{ can be expressed 
in terms of E\\, E 22 , vj 2 and Gn- 


Now the total deflection response can be evaluated as follows: 

+( £ ( SwL > Bbf ) b; +0*1 !dq{ ) ?; ) (2.76) 

1 

The variance of the response can be expressed by: 

Var( )=E[(S( do’L +0>*1 ) 2 ] (2.77) 

1 

The SD is the square root of the variance. The above approach is general in 
nature, and can be applicable to different forced vibration problems involving random 
system parameters and loading. 


25 



CHAPTER 3 


RESULTS AND DISCUSSION 

3.1 INTRODUCTION 

In this chapter numerical results have been presented and the behaviour of 
composite flat plates with random material properties and random loading has been 
described. This chapter is divided into three sections. First section describes the nonlinear 
behaviour of static deflection, the second section describes the nonlinear free vibration 
whereas the third section describes the nonlinear forced vibration behaviour. All edges 
simply supported boundary condition is used in all the three section. 

Closed form expressions have been developed for the mean and variance of 
response for simply supported plates in the previous chapter. These expressions are used 
to generate the numerical results. Results are obtained for symmetric and antisymmetric 
cross' ply laminates. All lamina are assumed to have the same thickness and the material 
properties are orthotropic along the material axes. The validation of the developed 
expressions is checked with available results in literature. The effects of material property 
dispersion along with variations in aspect ratio, thickness ratio and oscillation amplitude 
on the frequency statistics have been studied. 


3.2 STATIC DEFLECTION 
3.2.1 Validation 

The validation of the present formulation is sought by comparison of results with 
reported literature. The comparison has to be made with the linear formulation only, as 
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nonlinear formulation and results are not available for laminated composite plate with 
random material properties and loading. Based on equations (2.32), (2.33) and (2.35) 
deflection variance for linear strain-displacement relation obtained by taking L-f=Lz =0, 
with all basic material properties varying simultaneously are compared with the result by 
Salim et al. [10]. Table 3.1 represents the comparison of SD of center point deflection for 
four layered [0°/90 o /90 o /0°] square symmetric cross ply with thickness ratio b/h= =100. The 
mean value of the material properties and loading used are [10]: 

i?i i=l 81 GPa, £ 22 =10.3 GPa, Gi 2 =7.17GPa, v 12 =0.28and 20=1.65 x 10 3 Pa 
where E\\ and £ 22 denote longitudinal and transverse elastic modulus, respectively, Gi 2 is 
the inplane shear modulus and v )2 denotes the Poisson’s ratio. The plate edges are simply 
supported and the mid point deflection is non-dimensionlised with the plate thickness. A 
reasonable good agreement between the two is observed. The results from the present 
study with nonlinear formulation are also placed in the table for comparison. The present 
approach is expected to give more accurate results as it accounts for non-linearity in the 
system equation. Figure 3.1 presents the comparison of the linear response with nonlinear 
cases for different loadings. It can be seen that when loading is increased, the linear 
approach loses accuracy as compared to the present approach. 

3.2.2 Second Order Deflection Statistics 

The second order statistics of transverse deflection has been obtained for 
symmetric and anti-symmetric cross-ply rectangular composite plate with all edges 
simply supported. The laminates are subjected to a transverse sinusoidal load. For 
graphite/epoxy composite plate the mean values of the variables considered as random 
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for the analysis are - the material properties E U =25E 22 , G u = 0.5 E 22 , v u =0.25 and the 
external loading Q 0 = 2.0 x 10 3 Pa 

The influence of scattering in the material properties on the static deflection has 
been obtained by allowing the coefficient of variation to change from 0 to 20 % for 
laminated cross ply plates. The plates have aspect ratios b/a= 1, 2, 3 and 4 and thickness 
ratio 100 with stacking sequences of [0°/90°/9070°] and [0790°/0°/90°]. 

3.2.2.1 Mean Deflection - Table 3.2 presents the non-dimensionalised mean central 
deflection for the four-layered symmetric laminate with different aspect ratios. Some 
available results by Singh et al. [8] are also placed in the table for comparison. The 
reference uses Rayleigh-Ritz method for the solution. The present approach gives an 
exact solution for the mean analysis. A reasonable good agreement between the two is 
observed. 

Table 3.3 presents the non-dimensionalised mean value of the central deflection 
for different plate thickness with simply supported edges. The mean deflection, as 
expected, decreases with increase in the plate thickness. 

3.2.2.2 Variance of deflection - The variation of non-dimensionalised deflection with SD 
of all the basic inputs changing simultaneously for three cases-square symmetric, square 
anti-symmetric and rectangular symmetric four layered cross-ply with a/b=2 is 
represented in Figure 3.2. It is seen that the change in deflection for square symmetric 
cross-ply is greater than the other two cases. Figures 3.3 (a)-(d) present deflection 
sensitivity to dispersion in only one basic material variable at a time for bfh= 100. The 
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variations of deflection are most affected by changes in E \\ . It is observed that the 
deflection is least affected with dispersion in v\i. Figure 3.4 shows deflection sensitivity 
to dispersion in load only. It is seen that variation of load has a dominant effect on the 
deflection scattering as compared to E 22 , Gn and v\ 2 . 

Change in SD of deflection due to dispersion in all the basic random inputs 
changing simultaneously for symmetric cross-ply with different aspect ratio a/b= 1, 2, 3 
and 4 are presented in Figure 3.5. The effect on the scattering of deflection due to all 
basic random inputs decreases with increase in the aspect ratio. 

Figure 3.6 shows midpoint deflection sensitivity to dispersion in all the basic 
random inputs changing simultaneously for different plate thickness ratios (b/h)=100, 50 
and 33.33. The effect on the scattering of deflection increases with increase in thickness 
whereas for the linear case it is independent of variations in thickness. 

Figure 3.7(a)-(b) show variation in deflection SD with dispersion in all the basic 
random inputs changing simultaneously for square anti-symmetric and rectangular 
symmetric plates with different thickness ratios of 100, 50 and 33.33. Figure 3.8 (a)-(c) 
show deflection sensitivity to dispersion in all the basic random inputs changing 
simultaneously with different loading for square symmetric, square antisymmetric and 
rectangular symmetric cross ply. A comparison with linear analysis with same boundary 
condition and loading is also presented. It is observed that the effect on scattering of 
deflection decreases with increase in loading whereas for linear case it is insensitive to 
change in mean loading. 
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Table 3.1: Comparison of SD/h of center point deflection for [0 O /90°/90 < 70°] laminate 
with a/b=l and thickness ratio b/h=\00 


S D/mean of 
material properties 

0.05 

0.10 

0.15 

(SD/h, wo) 




Salim [6] (Linear ) 

0.00360 

0.00715 

0.01080 

Present Study 
(Linear) 

0.00371 

0.00724 

0.01091 

Present Study 
(Non-Linear) 

0.00328 

0.00690 

0.01003 


Table 3.2: Nondimensionalised mean central deflection 

laminate with different aspect ratios and Mz=100 

(wm/A) for [0°/90 O /90°/0 O ] 



Aspect Ratio (a/b) 



a/b= 1 

a/b = 2 afb = 3 

a/b=A a/b=5 

Present work 

0.518 

0.8124 0.8435 

0.8503 0.8527 

Ref [13] 

0.521 

0.8128 
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Table 3.3: Nondimensionalised mean value of the central deflection (wn\/h) for different 
plate thickness 



Thickness ratio (b/h) 




100 

50 

33.33 

25 

(a) Stacking 

sequence: [0°/9079070 u ] : 




a/b=l 

0.5180 

0.1069 

0.0319 

0.0135 

a/b= 2 

0.8124 

0.4077 

0.1493 

0.0638 

(b) Stacking 

sequence: [o7907o790°] : 




a/b= 1 

0.5431 

0.1250 

0.0374 

0.0158 
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sd/mean, 



Figure 3.1- Comparison for coefficients of variation of deflection of [0 o /9O°/90 o /0°] 
laminate having a/b - 1 and b/h=\00 with all material properties changing simultaneously. 
Casel - g 0 =1.65 x 10 3 Pa, Case2 - Q 0 =2.65 x 10 3 Pa, Case3 - Q 0 =3.65 x 10 3 Pa, Case4 - 
<2o=4.65 x 10 3 Pa 
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sd/mean, w. 



Figure 3.2- Influence of SD of all basic random inputs changing simultaneously on 
coefficient of variation of deflection for different cross-ply with b/h= 100. 
CaseI-[0°/90 o /90 o /O°] laminate with a/b=l, Casen-[0 o /9O o /0°/90 o ] laminate with a/b= 1, 
Casein- [O 0 /9O°/9O°/O 0 ] laminate with a/b= 2.0 


33 




caselll 



Figure 3.3- Influence of SD of basic material properties on coefficient of variation of 
deflection of different cross ply with b/h= 100, (a) only En varying, (b) only E 22 varying, 
(c) only G 12 varying, (d) only v 12 varying. (Three cases same as in Figure 3.2) 








sd/mean, 



Figure 3.4- Influence of SD of load l q k on coefficient of variation of deflection for 
different cross ply with Mz=100. (Three cases same as in Figure 3.2) 
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sd/mean, 



Figure 3.6- Influence of SD of all basic random inputs changing simultaneously on 
coefficient of variation of deflection for [0°/90 o /90°/(r] laminate having different plate 
thickness with afb = 1 
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Figure 3.7- Influence of SD of all basic random inputs changing simultaneously on 
coefficient of variation of deflection of different cross-ply with different plate thickness 
for (a) [0 0 /90°/0 0 /90°] laminate with a/b= 1, (b) [0 0 /9O°/9O°/O 0 ] laminate with a/b = 2 
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Figure 3.8- Influence of SD of all basic random inputs changing simultaneously on 
coefficient of variation of deflection of different cross-ply having different loading on the 
plate with b/h=\00. (a) [0°/90 0 /90°/0 0 ] laminate with a/b=l, (b) [0°/90°/0 0 /90 o ] laminate 
with a/&=l, (c) [O°/9O 0 /90 0 /O 0 ] with a/b= 2; loading 1- Q 0 = 2.0 x 10 3 Pa, loading 2- Q 0 = 4.0 
x 10 3 Pa. loading 3- 6.0 x 10 3 Pa, loading4 - 8.0 x 10 3 Pa. 






3.3 FREE VIBRATION 


3.3.1 Validation 

The validation of the present formulation is sought by comparison of results with 
reported literature. However, nonlinear formulation is not available in literature for 
laminated composite plate with random material properties hence comparison has been 
made with the linear formulation only with the linear form as a special case of the present 
nonlinear formulation. Based on equations (2.54) and (2.56) frequency variance for linear 
strain-displacement relation obtained by taking L 3 - 0, with all basic material properties 
varying simultaneously are compared with the result by Singh et al. [25], Table 3.4 
represents the comparison of SD of frequency for four layered [079079 070°] square 
symmetric cross ply with thickness ratio b/h=l0. The mean values of the material 
properties used are [25]: 

£■22=10.3 GPa, Eu= 25E 2 2, Gu^O-SEn, and vi 2 = 0.25 
The plate edges are simply supported and the SD of frequency is non-dimensionalised 
with the mean frequency. A reasonable good agreement between the two is observed. The 
results from the present study with nonlinear formulation are also placed in the table for 
comparison. The effect of nonlinearity is apparent with different value of the amplitude 
Wm<u* 

Table 3.5 presents a comparison of the non-dimensionalised mean frequency for 
the four-layered symmetric laminate with aspect ratio a/b= 2 having different amplitudes 
with limited available results by Singh et al. [14]. The material properties used for the 
analysis are [14]: 

E\ \ =40E 22 , Gi2=0.5£ 2 2, and v 12 =0.25. 
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The reference uses direct numerical integration method for the analysis of mean 
frequency whereas the present approach gives an exact solution. A reasonably good 
agreement between the two is observed. 

3.3.2 Second order frequency statistics 

The material used for the graphite/ epoxy composite plate is the same as employed 
for generating Table 3.5. All the four material properties are considered as random for the 
analysis. 

3.3.2. 1 Mean frequency 

Table 3.6 presents the non-dimensionalised mean nonlinear frequency for 
different plate thickness and amplitudes with simply supported edges. Influence of the 
scattering in the material properties on the mean frequency has been obtained by allowing 
the coefficient of variation to change from 0 to 20 % for laminated cross ply plates. The 
plates have aspect ratios b/a = 1 and 2 and thickness ratios b/h = 100, 50 and 33.33 with 
stacking sequences of [0 o /90°/90 o /0°] and [0°/90 o /0° 90°]. The mean nondimensional 
frequency decreases with increase in the plate thickness and increases with the oscillation 
amplitude. The anti-symmetric laminate has higher mean frequency compared to the 
symmetric laminate. 

3.3.2.2 Variance of frequency 

The variations of non-dimensionalised frequency with dispersion in all the basic 
material properties changing simultaneously are presented in Figures 3.9(a) - (c). Three 
cases square symmetric, square anti-symmetric and rectangular {aJb— 2) symmetric four 
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layered cross ply with fofr-100 have been examined. It is seen that as normalised 
amplitude (w’=100 w max /b) increases the frequency sensitivity increases for all cases. The 
sensitivity of the square anti-symmetric cross ply is greater than the symmetric case. The 
frequency dispersion also increases with increase in the aspect ratio. 

Figures 3.10 (a) - (d) present the frequency sensitivity to dispersion with only one 
basic variable random at a time with different amplitudes for square symmetric cross ply 
for b/h=\00. The frequency variations are most affected by change in E n and least 
affected by dispersion in vn. It is also seen that sensitivity of the plate frequency to 
oscillation amplitude increases with increased variations in E u and decreases with En, v t2 
and Gn- Figures 3.11 (a) - (d) and Figures 3.12 (a) - (d) show the frequency sensitivity to 
dispersion for only one basic variable random at a time with different amplitudes for 
square antisymmetric and rectangular symmetric cross ply respectively for b/h= 100. 

Frequency sensitivity to dispersion in all the basic random inputs changing 
simultaneously for symmetric cross-ply with aspect ratio a/b= 1, 2, 3 and 4 are presented 
in Figure 3.13(a)-(d). The effect on the scattering of frequency increases with increase in 
the aspect ratio from 1 to 2, decreases slightly from 2 to 3 and then remains almost 
constant for aspect ratios greater than 3. 

Figure 3.14(a)-(c) shows changes in frequency SD due to dispersion in all the 
basic random inputs changing simultaneously for square symmetric cross-ply with 
different amplitude and plate thickness ratios (b/h)=100, 50 and 33.33. The effect on 
frequency scatter shows an initial decrement with increase in thickness ratio. This effect 
is less pronounced for thicker plates. For the linear case, frequency is independent of 
variations in both thickness ratios and oscillation amplitudes. 
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Figures 3.15(a)-(c) and 3.16(a)-(c) show frequency sensitivity to dispersion in all 
the basic random inputs changing simultaneously for square antisymmetric and 
rectangular symmetric cross-ply with different amplitude and plate thickness ratios 
(b/h)— 100, 50 and 33.33. The effects on frequency scatter for both square anti-symmetric 
and rectangular symmetric laminate due to thickness ratio has the same nature as that of 
the symmetric cross-ply. 

Influence of amplitude (w mar /h) on frequency coefficient of variation with a 
constant SD/mean for all material properties changing simultaneously for symmetric 
cross-ply with b/h=\00 is shown in Figure 3.17. It is found that the frequency scatter 
increases nonlinearly with increase in amplitude and also increases with increase in 
variation of material properties. 

3.4 FORCED VIBRATION 
3.4.1 Validation 

The validation of the present formulation is sought by comparison of results with 
reported literature. However, nonlinear formulation is not available in literature for 
laminated composite plate with random material properties and random loading hence 
comparison has been made with the mean analysis only. 

Table 3.7 presents a comparison of the non-dimensionalised mean frequency for 
the square two-layered cross-ply laminate with different amplitude with limited available 
results by Singh et al. [19]. The material properties used for the analysis are [19]: 

E\ i =40£’ 2 2, Gi2=0.5£ 2 2, and v 12 =0.25, and excitation Q 0 =0.2 L\ 
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Table 3.4: Comparison of SD( ® ! )/mean( ) f„ r [0°/90°/9oV] laminate with a/b= 1 and 
thickness ratio b/h= 10 


SD/mean of 
material properties 

0.05 

0.10 

0.15 

0.20 

(SD( 03 2 )/mean ( a 2 ) 





Singh [6] (Linear) 

0.045 

0.091 

0.138 

0.180 

Present Study 
(Linear) 

0.044 

0.088 

0.133 

0.177 

Present Study (w’=0.3) 
(Non-Linear) 

0.048 

0.097 

0.146 

0.194 

(w’=0.6) 

0.051 

0.103 

0.154 

0.207 

II 

O 

kO 

0.054 

0.109 

0.162 

0.221 


Table 3 . 5 : Comparison of the nondimensionalised mean nonlinear frequency ( & ni/ <z> i ) for 

[0 O /90 O /90°/0°] laminate with aspect ratio a/b= 2 and Mi=100 


>i;/nax=0.3 

Wmax= 0.6 

Wmax= 0.9 

Present work 

1.18 

1.61 

2.16 

Ref [13] 

1.22 

1.63 

2.18 
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Table 3.6: Nondimensionalised mean nonlinear frequency (®ni/©i) for different plate 
thickness and amplitudes 


Thickness ratio ( b/h ) 


100 50 33.33 


(a) Stacking sequence: 

[0 /9079070 u ]: 



a/b= 1 

(w 5 =0.3) 

1.082 

1.021 

1.009 


(w’=0.6) 

1.300 

1.082 

1.037 


(w’=0.9) 

1.602 

1.178 

1.082 

a/b - 2 

(w’=0.3) 

1.204 

1.054 

1.024 


(w’=0.6) 

1.680 

1.204 

1.095 


(w’=0.9) 

2.265 

1.421 

1.204 

(b) Stacking sequence: 

[0°/90°/0 O /90 0 ]: 



a/b= 1 

(w’=0.3) 

1.097 

1.025 

1.011 


(w’=0.6) 

1.351 

1.097 

1.044 


(w’=0.9) 

1.695 

1.209 

1.097 
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Figure 3.9- Influence of SD of all basic random inputs changing simultaneously on 
coefficient of variation of frequency of different cross-ply with b/h= 100. 

(a) [0 o /90°/9O o /0°] laminate with a/b=\, (b) [O o /9O°/O 0 /90 0 ] laminate with a/b= 1, (c) 
[O o /9O o /90 o /O°] laminate with a/b-2.0 
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Figure 3.13- Influence of SD of all basic random inputs changing simultaneously on 
coefficient of variation of frequency for [0°/90 o /90 o /0 0 ] laminate having different aspect 
ratio of plate, (a) a/b=\, (b) a/b- 2, (c) a/6=3 and (d) a/6=4 
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coefficient of variation of frequency for [0°/90°/0°/90°] laminate having different plate 
thickness with a/b= 1. (a) b/h-lOQ, (b) b/h=50 and (c) b/h— 33.33 
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Figure 3.17- Influence of amplitude 
[0790°/90°/0°] laminate with all materia! 


5 - 



The reference uses direct numerical integration method for the analysis of mean 
frequency whereas the present approach gives an exact solution. A reasonably good 
agreement between the two is observed. 

3.3.2 Second order amplitude statistics 

The material used for the graphite/epoxy composite plate is the same as employed 
for generating Table 3.7. All the four material properties and excitation are considered as 
random for the analysis. 

3.3.2.1 Mean amplitude 

Table 3.8 presents the non-dimensionalised mean nonlinear frequency for 
different plate thickness and amplitudes with simply supported edges. Influence of the 
scattering in the material properties on the mean amplitude has been obtained by allowing 
the coefficient of variation to change from 0 to 20 % for laminated cross ply plates. The 
plates have aspect ratios b/a - 1 and thickness ratios bfh = 100, 50 and 33.33 with stacking 
sequences of [0 o /90°/90 o /0°] and [0 o /90°/0°/90 o ]. The mean amplitude decreases with 
increase in the plate thickness and increases with the oscillation amplitude. The anti- 
symmetric laminate has higher mean amplitude compared to the symmetric laminate. 

3.3.2.2 Variance of amplitude 

The variations of non-dimensionalised amplitude with dispersion in all the basic 
inputs changing simultaneously are presented in Figures 3.18 and 3.19. Two cases square 
symmetric and anti-symmetric four layered cross ply with b/h=\Q0 have been examined. 
It is seen that as excitation frequency increases the amplitude sensitivity increases for all 
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cases. The sensitivity of the square anti -symmetric cross ply is greater than the symmetric 
case. Here, the graphs are presented for excitation frequencies ©1=31.62, © 2 =47.43 and 
©3=63.34 rad/sec. 

Figures 3.20 (a) - (d) and 3.21 (a) - (d) present the influence of only one basic 
material property random at a time on displacement response dispersion for square 
symmetric and antisymmetric cross ply with M=100. The amplitude variations are most 
affected by change in E n and least affected by dispersion in v i2 . It is also seen that 
sensitivity of the oscillation amplitude to plate frequency increases with increased 
variations in En, En, Gn and vi 2 . 

Figures 3.22 and 3.23 show amplitude sensitivity to dispersion in excitation for 
square symmetric and antisymmetric cross-ply respectively for b/h=l00. It is seen that 
variation of excitation has a dominant effect on the deflection scattering as compared to 
En, Gn and v n . 

Influence of frequency on coefficient of variation of amplitude with a constant 
SD/mean for all basic inputs changing simultaneously for square symmetric cross-ply 
with b/h-\00 is shown in Figure 3.24. It is found that the amplitude scatter increases 
nonlinearly with increase in frequency and also increases with increase in variation of 
basic inputs. 

The interrelation between (ww/h) on coefficient of variation of amplitude with a 
constant SD/mean for all basic inputs changing simultaneously for square symmetric 
cross-ply with b/h = 100 is shown in Figure 3.25. This also shows that the amplitude 
scatter increases nonlinearly with increase in amplitude and also increases with increase 
in variation of basic inputs. 
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Table 3.7 : Comparison of the nondimensionalised mean nonlinear frequency ( co n |/ co |) for 
[0790°] laminate with different amplitude and aspect ration a/b=\ 


WmJh 

Present work 

Ref[19 ] 

+0.4 

0.8943 

0.8898 

-0.4 

1.3408 

1.3395 

+0.6 

1.1595 

1.1486 

-0.6 

1.4167 

1.4109 

+1.0 

1.6428 

1.6116 

-1.0 

1.7581 

1.7334 


Table 3.8: Nondimensionalised mean nonlinear frequency (coj a>\) for different square 
plates with different amplitudes 


(&) n \/ CO |) 


w max /h 

Stacking sequence 
[079079070°] 

Stacking sequence 
[079070790°] 

+0.4 

0.8984 

0.9311 

-0.4 

1.3435 

1.3654 

+0.6 

1.1667 

1.2233 

-0.6 

1.4224 

1.4689 

+1.0 

1.6568 

1.7676 

-1.0 

1.7718 

1.8749 
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Figure 3.20- Influence of SD of basic material properties on coefficient of variation of 
amplitude for [0°/90 0 /0°/90 0 ] laminate with a/b= 1.0 and b/h= 100. (a) only Eu varying, (b) 
only E 22 varying, (c) only Gn varying, (d) only v ]2 varying 
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Figure 3.21- Influence of SD of basic material properties on coefficient of variation of 
amplitude for [0 0 /90°/90°/0 0 ] laminate with a/b= 2.0 and b/h=l00, (a) only En varying, (b) 
only E 22 varying, (c) only G 12 varying, (d) only v \ 2 varying 
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Figure 3.22- Influence of SD of excitation ‘ q *’ on coefficient of variation of amplitude 
for [O o /90°/90 o /O°] laminate with a/b = 1 .Oand b/h= 100 













CHAPTER 4 


CONCLUSIONS 

An approach has been presented in Chapter 2 to obtain the second order statistics 
for static deflection, free vibration and forced vibration response of laminated composite 
plates with random material properties under random loading for nonlinear strain- 
displacement relations. The system equation for the plate has been solved using termwise 
series integration and a first order perturbation technique is used to obtain the second 
order statistics. 

Results obtained using the present approach have been discussed in Chapter 3 for 
two stacking sequences of cross-ply laminates. The following conclusions are noted from 
the results for the graphite-epoxy laminated plate having all edges simply supported with 
movable edges. 

Static deflection 

(i) The dispersions in the static deflection show linear variation with SD of the material 
properties and the external loading in the range studied. One of the main reason for 
this may be the FOPT employed in developing the analysis. 

(ii) The influence of SD of static deflection shows different sensitivity to different input 
variables. The sensitivity changes with the laminate construction, a/h ratios and the 
material. However, any recognizable pattern is not apparent in all cases. 

(iii) Increase in thickness ratio of the plate results in increase in deflection scatter whereas 
in case of linear analysis plate thickness ratio has no effect. 

(iv) Increase in the aspect ratio of the plate results in decrease in deflection scatter. 
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(vi) Variation of E\\ and load Qo has dominant effect on the scattering of deflection as 
compared to E22, G n and v n . 

(vii) Variation in the external load causes larger variations in the deflection response 
compared to the effect of material properties. 

(viii) Increase in mean value of load results in decrease in deflection scatter whereas it is 
predicted to be constant by linear analysis. 

Free vibration 

(i) The dispersions in the frequency response show linear variation with SD of the 
material properties in the range studied whereas nonlinear variation in frequency has 
been seen with variation in the oscillation amplitude. 

(ii) Increase in thickness ratio and oscillation amplitude of the plate results in increase in 
frequency scatter with all material properties changing simultaneously. Linear 
analysis shows no dependence of frequency on plate thickness ratio and oscillation 
amplitude. 

(iii) Variation of E\\ has dominant effect on the scattering of frequency as compared to 
E22, G\2 and V\2- 


Forced vibration 

(i) The dispersions in the response amplitude show linear variation with SD of the 
material properties in the range studied whereas nonlinear variation in amplitude has 
been seen with variation in the oscillation amplitude and frequency. 
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(ii) Variation of E u has dominant effect on the scattering of amplitude as compared to 
E22, Gn and V \ 2 and the excitation Qo. 

(iii) Variation in the external load causes larger variations in the deflection response 
compared to the effect of material properties. This is also observed in the static case. 

SUGGESTIONS FOR FUTURE WORK 

The present study gives a first attempt at nonlinear analysis for laminated composite 

plates with random material properties subjected to random loading. The following 

extensions can be applied to the present problem. 

(i) The problem can be solved by using higher order perturbation technique by taking 
higher order terms in Taylor series expansion of material properties and loading to get 
higher dispersion. 

(ii) Effect of correlation between the material properties as well as between material 
properties and loading on the response should be studied. 

(iii) Material properties can be assumed to vary with respect to space and time to 
incorporate the aging effect as well as point to point variations in composites. 

(iv) Random loading can be assumed to vary with respect to space also. 

(v) Geometric parameters and boundary condition can also be modeled as random. 

(vi) Effect of FSDT and HSDT on the response can be examined to extend the 
formulation to thick plates. 

(vii) Effect of hygrothermal behaviour on the response can also be examined. 

(viii) The formulation can be extended to curved panels. 
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